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Abstract 

The theory for Bose-Einstein correlations in case of neutral gauge bosons in pp colli- 
sions at high energies is presented. Based on quantum field theory at finite tempera- 
ture the two-particle Bose-Einstein correlations of neutral gauge bosons are carried out 
for the first time. As a result, the important parameters of the correlation functions 
can be obtained for the Z°Z° pairs. The correlations of two bosons in 4-momentum 
space presented in this paper offer useful and instructive complimentary viewpoints 
to theoretical and experimental works in multiparticle femtoscopy and interferometry 
measurements at hadron colliders. 



1 Introduction 



An investigation of the space-time extension or even squeezing of particle sources 
via the multiparticle quantum-statistics correlation in high energy interactions is still 
attract the attention of physical society in both experiment and theory. Over the 
past few decades, a considerable number of successful studies have been done in this 
direction [1] . It is well understood that the studies of correlations between produced 
particles, the effects of coherence and chaoticity, an estimation of particle emitting 
source size play an important role in high energy physics. 

By studying the Bose-Einstein correlations (BEC) of identical particles (we mean 
like-sign charge particles and the neutral charge ones), it is possible experimentally 
to determine the time scale and spatial region over which particles do not have the 
interactions. Such a surface is called as decoupling one. In fact, for an evolving system 
such as pp collisions, it is not really a surface, since at each time there is a spread 
out surface due to fluctuations in the last interactions, and the shape of this surface 
evolve even in time. The particle source is not approximately constant because of 
energy-momentum conservation constraint. 

More than half a century ago Hanbury-Brown and Twiss [2] used BEC between 
photons to measure the size of distant stars. In the works [3,4], the master equations 
for evolution of thermodynamic system that can be created at the final state of a high 
multiplicity process were established. The equations have the form of the field oper- 
ator evolution equation (Langevin-like [5]) and allows one to gain the basic features 
of the emitting source space-time structure. In particular, it has been conjectured 
and further confirmed that the size of the emitting source through BEC is strongly 
affected by non-classical off-shell effect. 

The shapes of BEC function were experimentally established in the LEP experi- 
ments ALEPH [6], DELPHI [7] and OPAL [8], and ZEUS Collaboration at HERA [9], 
which also indicated a dependence of the measured correlation radius on the hadron 
(-7T, K) mass. The results for 7r ± 7r ± and ir ± ir T correlations with pp collisions at y/s 
= 1.8 TeV were published by E735 Collaboration in [10]. 

The correlations between heavy particles (e.g., neutral gauge bosons Z a Z°) of 
Bose-Einstein type have not been carried out previously at hadron colliders. Such a 
study can be addressed to the Large Hadron Collider (LHC) which provide proton- 
proton interactions at \fs = 14 TeV centre-of-mass system (c.m.s.) energy. 

In this work, we make an attempt to demonstrate that the problem of properties 
of the genuine interactions can be explored using experimental data which can be 
collected by ATLAS and CMS Collaborations at the LHC. These data can be ana- 
lyzed through the compared measures of some inclusive distributions and final state 
correlations. 

One of the aims of this paper is to carry out the proposal for the experimental 
measurements of the Z°Z° pair correlations. 

This exploration will be theoretically supported by the quantum field theory at 
finite temperature (QFTp) model approach [3]. It is known that the effective tem- 
perature of the vacuum or the ground state or even the thermalized state of particles 
distorted by external forces is occurring in models quantized in external fields. One of 
the main parameters of the model is the temperature of the particle source under the 
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random source operator influence. The main channels are the di-lepton production 
pp — ► Z°Z° — > 2e~2e + , 2/i~2/x + , e~e + ^i~/i + in collisions. 

An efficient selection of leptons needs to be performed according to the following 
criteria. First, all leptons were required to lie in the pseudorapidity range covered 
by, e.g., the CMS muon system that is, \r)\ < 2.4. Second, the leptons were required 
to be unlikely charged in pairs. Note that the acceptances of another multipurpose 
detector ATLAS in the azimuthal angle and pseudorapidity are close to the respective 
parameters of CMS. 

The dilepton channel is especially promising from the experimental point of view, 
since it is expected that the experimental facilities related for LHC (CMS and ATLAS 
detectors) will make it possible to record muons of energy in the TeV range with 
a resolution of about a few percent and an efficiency close to 100 %. Moreover, 
this channel is characterized by a maximum signal-to-background ratio in the energy 
region being considered. 



2 BEC in case of two particles 

A pair of identical bosons with the mass m produced incoherently (in ideal nondis- 
turbed, noninteracting cases) from an extended source will have an enhanced proba- 
bility Ciipx.p-i) — Ni2(pi,P2)/[Ni(pi)- N2(j>2)] to be measured in terms of differential 
cross section a, where 

to be found close in 4-momcntum space 3?4 when detected simultaneously, as compared 
to if they are detected separately with 

1 da d 3 pi 



N ^=aWk' * = W^-' £ -V P « +m ' i=1 ' 2 - (2) 
On the other hand, the following relation can be used to retrieve the BEC function 

where N(Q) in general case refer to the numbers for neutral gauge bosons (eg., Z°Z°) 
with 



Q = y/-(pi - p 2 ) M • (pi - V2Y = V / M 2 -4m 2 . (4) 
In definitions © and (JD), N re f is the number of particle pairs without BEC and 



p lli — (u)i,pi) are four-momenta of produced bosons (i — 1, 2); M = y (pi + P2)f l is 

the invariant mass of the pair of bosons. 

An essential problem in extracting the correlation is the estimate of the reference 
distribution N re f(Q) in Eq. If there are other correlations beside the Bose- 

Einstcin effect, the distribution N re f(Q) should be replaced by a reference distribution 
corresponding to the two-particle distribution in a geometry without BEC. Hence, the 
expression ([3]) represents the ratio between the number of Z°Z° pairs N(Q) in the 
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real world and the reference sample N re f (Q) in the imaginary world. Note that the 
reference sample can not be directly observed in an experiment. Different methods 
are usually applied for the construction of reference samples [1], however all of them 
have strong restrictions. One of the preferable methods is to construct the reference 
samples directly from data. For our aim for reference sample N re f(Q), it is suitable 
to use the pairs Z°Z° from different (mixed) events. 

It is commonly assumed that the maximum of two-particle BEC function C2(Q) is 
2 for pi — P2 if no any distortion and final state interactions are taking into account. 

There are experimental difficulties in a determination of Z°Z° correlations, which 
are associated with acceptance limitations and limited statistics in the Z°Z° sample. 

In general, the shape of the BEC function C^Q) is model dependent. The most 
simple form of Goldhabcr-likc parameterization for C2(Q) [11] has been used for data 
fitting: 

C 2 (Q) = Co-(l + \e~ Q2R2 )-(l + eQ), (5) 

where Co is the normalization factor, A is the so-called the chaoticity strength factor, 
meaning A — 1 for fully incoherent and A = for fully coherent sources; the parame- 
ter R is interpreted as a radius of the particle source, often called as the "correlation 
radius" , and assumed to be spherical in this parameterization. The linear term in 
(JS|) is supposed to be account within the long-range correlations outside the region of 
BEC. Note that distribution of bosons can be either far from isotropic, usually con- 
centrated in some directions or almost isotropic, and what is important that in both 
cases the particles are under the random chaotic interactions caused by other fields 
in the thermal bath. In the parameterization ([5]) all of these problems are embedded 
in the random chaoticity parameter A. To advocate the formula ([5]) it is assumed: 

a. incoherent average over particle source where A serve to account for: 

- partial coherence, 

- long-lived resonances associated with multiple distinguishable sources, 

- Z°Z° purity; 

b. spherical Gaussian density of particle emission cell (with radius R); 

c. static source which means no time (energy) dependence. 

In order to save the quantum pattern of particle production process and to avoid 
the static and undistorted character of particle emitter source we also suggest to use 
the C2{Q) function within QFTp accompanying by quantum evolution approach in 
the form: 



r I - 2\/aa' P- 1 

c- 2 (Q) = OT-[i + (1 + a)(1 + a/) "(Q)+ (1 + a ) (1 + a0 ^(Q)]-F(Q,A X ), (6) 

where £(N) depends on the multiplicity N as 

m - or, 

The function F(Q, Ax) that expresses the correlation magnitude as a function of Q 
and two-particle relative distance Ax is a consequence of the Bogolyubov's principle 
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of correlations weakening at large distances [12] 



The function |jHJ) is normalized as F(Q,Ax = oo) = 1, and rt is the measure of 
correlations weakening where rj — ► as Ax — ► oo; f(Q,Ax) is the two-particle 
distribution function with Ax, while f(pi) are one-particle probability functions with 
£=1,2. 

The important parameter a (as well as a') in (JS]) summarizes our knowledge of 
other than space-time characteristics of the particle emitting source, and plays the 
role of a coherence parameter (see [4] for details). 

The f2(g) in ([5]) has the following structure in momentum space 

O(Q) = fi(Q) • 7 (n), (9) 

where 

fi(Q) = exp(-A pSR ) = exp - p 2 )^ ( Pl - j*)"] (10) 

is the smearing smooth dimensionless generalized function, is the (nonlocal) 
structure tensor of the space-time size (BEC formation domain), and it defines the 
spherically-like domain of emitted (produced) bosons. 

To clarify with 7(n) in @ let us emphasize that most of experiments dealing with 
elementary particles at high energies are of an inclusive as one measures quantum 
effect of BEC on limited samples of particles produced only. The unobserved part 
of the rest particle system acts then as a kind of thermal (heat) bath influencing 
measured samples of data (observables). Actually, the temperature T being the most 
important parameter describing the influence of such a thermal bath is occurred in 
this model. 

The function 7(71) reflects the quantum thermal features of BEC pattern and is 
defined as 

n 2 (ui) 1 u> + oj' 

where n(u, (3) is the mean value of quantum numbers for Bose-Einstein statistics parti- 
cles with the energy lu and the chemical potential fi in the thermal bath with statistical 
equilibrium at the temperature T = 1//3. The following condition Ylf n f( w i 0) = N 
is evident, where the discrete index / reflects the one-particle state /. 

In terms of time-like Rq, longitudinal Rl and transverse Rt components of the 
space-time size R^ the distribution A p jf looks like: 

ApSR - A pR = (Ap°) 2 R 2 Q + {Ap L fRl + (Ap T ) 2 R T . (12) 

Seeking for simplicity one has (Rl = Rt = R) 

\r - (p§ - P^fRl + (Pi - P2?R 2 (13) 

for identical bosons. 
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Hence, we have introduced a new parameter R^, a 4- vector, which defines the 
region of nonvanishing particle density with the space-time extension of the particle 
emission source. Expression (fTUf must be understood in the sense that Q(Q) is a 
function that in the limit R — > oo, strictly becomes a 5-function. For practical using 
with ignoring the energy- momentum dependence of a, and assuming that a 1 — a (a 
is related with 6*2(0) and N), we get the expression with ft(Q) ~ exp(— Q 2 R 2 ): 

C 2 {Q) ~ £(iV) {l + A„ ei0 (/3) e-Q 2R2 [l + A corr (/3) e +Q 2 « 2 /2] j ; (14) 

where the new intercept function becomes as A„ etu = 7(0;, /3)/(l + a) 2 , and the new 
coherence correction in the brackets of Eq. (I14[) carries an additional intercept func- 



tion X corr = 2 a/y/j(u), (3). In fact, since a 7^ a' (because u> ^ u>' and, therefore, 
the number of states identified here with the number of particles n{w) with given 
energy is also different), one can use the general precise form @ for Ci with details 
given by Eqs. (JH]) and pip and with a coherence function depending on the particle 
mass, the energy of final leptons produced in pairs within the decays of Z 0, s and such 
characteristics of the emission process as the temperature T and chemical potential 
/j. occurring in the definition of n(oS) in 

Since we did not follow special assumptions on the quantum operator level for 
C2{Q) from the initial stage, it may correspond to a physically real and observable 
effect at the LHC. This pattern may lead to a new squeezing state of correlation 
region. 



3 Stochastic field and Green's function 

Let us consider the stochastic field B^ix) = B^ g (x,r) that depends on the arbitrary 
random source s(x), and the fifth component r means the "stochastic time". The 
differential equation of an evolution of the field operator B^ M (x, r) = B^(x, r) in the 
system under the associated stochastic process is 



d T B fl (x,T)=0[B fi (x,T)}, (15) 
where 0[B^(x,t)] is the differential stochastic operator which has the form 

0[B,(x,r)} = -l ^Bfcll] + S,(x,r) (16) 
V SB^x^t) 

with a volume V being introduced by dimensional reason. The r.h.s. of Eq. (|16p 
is the so-called stochastic operator where J = f d 4 yL[B tl {y),d v B^ v {y)\ is the action 
defined by the Lagrangian density L; s^x^t) = s^{x,t) + n^P carries the random 
stochastic history where the memory dissipation forces and the heat bath effects are 
included into s^(x) — s^(x, r), the constant P emerges within the action of the 
stationary forces. Equation (|15| is nothing other but the evolution equation of the 
Langevin type applied already to stochastic processes on the quantum operator level 
in derivation of multiparticle Bose-Einstein correlations [3] . 
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For simplicity, we assume that s^(x) varies stochastically with the Gaussian cor- 
relation function 

(s fl (x)s u (y)) = const S^exp^z 2 /l 2 ch ) , (17) 

where z v — (x—y) v , and const is the strength of the noise described by the distribution 
function exp(— z 2 /l 2 h ) with l c h being the noise characteristic scale. Both const and 
l c h define the influence of the (Gaussian) noise on, e.g., correlations between particles 
that "feel" an action of an environment. Actually, Eq. (fT5|) can be transferred to the 
standard field equation of motion (in Euclidean space) 

(18) 



V SB„(x) 



with the source s M (x) if both and s M do not depend on "stochastic time" r. In 
classical theory, the random process given by s^^x) is nothing other but the white 
(Gaussian) noise. 

In this paper, we focus on the role of particle masses and energies, effects of 
coherence and distortion, and the heat bath influences which arc rather important to 
describe the correlations between particles. To solve this problem, especially to derive 
the memory term in evolution equation one can use the general properties of QFTp. 
The model is defined by the following generating functional in four-dimensional space- 
time 



Z = 

where 



J DB^cxp -i J d A xL{B^,B^) 



(19) 



■]b^B^ + l(m 2 + U)B^ (20) 



with B^y = d^By - ByB^. 

The direct calculations using the solution of Eq. (fl3)l with the Lagrangian density 
([20]) leads to the propagator of the field £? M (x, t) distorted by s fl (x,T). The trans- 
verse part of B )1 (x,t) will give the correct expression for the Euclidean vector field 
propagator at r — > oo. 

We are working with fields that correspond to a thermal field B fJi (x) with the 
standard definition of the Fourier transformed propagator F[G ^ v {p)\ 

F[G MV (p)] = G^(x-y) = Tr{T[B^x)B v (y)]p p } (21) 

with pp = e~@ H /Tre~@ H being the density matrix of a local system in equilibrium 
at temperature T under the Hamiltonian H 

« = j W<'°tf ^ <22) 

with the operators of annihilation b x (p) and creation 6 A+ (p) to be defined later. 
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The interaction of B^(x) with the external field is given by the potential U. The 
equation of motion is 

(V 2 + m 2 )-B M (x) = -J M (x), (23) 

where J^{x) — UB^(x) is the source density operator. A simple model like this allows 
one to investigate the origin of the unstable state of the thermalized equilibrium 
in a nonhomogeneous external field under the influence of source density operator 
Jfi(x) = UB^(x). For example, the source can be considered as 5-like generalized 
function, J^{x) = p, p(x,e)B ll (x) in which p(x, e) is a <5-like succession giving the 6- 
function as e — > (where p is some massive parameter) . This model is useful because 
the J-like potential U(x) provides the conditions for restricting the particle emission 
domain (or the deconfinement region). We suggest the following form: 

J M (x) = - J S ys(x) B M (x) + Jr^x), (24) 

where the source J^{x) is a sum of a regular systematic motion part J sys {x) and the 
random source Jr m (x). The equation of motion (f2"3"|) becomes 

[V 2 + to 2 - J sys }B^(x) = -J Rfl (x), (25) 

and the propagator satisfies the following equation (in the Fourier transformed form 
labeled by tilde): 

\?l ~ m2 + Jsys]G^(p^) = d^(p), (26) 

where 

1 d 2 



d ^ = {^ + ^8^dx-J 5{x) - (2?) 

As the standard point, the Green's function of the vector field can be obtained from 
the one of the scalar field acting by the relevant operator + m~ 2 d 2 /(dx^dx^). 
The solution of Eq. is 



B^x) = -J dyG^{x,y)J R Mi (28) 
where the Green's function obeys the Eq. (|26|) . 



4 Green's function and kernel operator 

Let us go to the thermal field operator B^ix) by means of the linear combination of 
the frequency parts B^(x) and -B 2+ (x) 

B^x) = Bl(x) + Bf (x) (29) 

with [B M (x), B v (y)] = iD^(x - y) and 

d 3 p 
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The following properties of polarization vectors are the standard ones: 

E^(^( P ) = -.v + ^. 

A=l 

We assume that the deviation from the asymptotic free state given by the operator 
a(p, t) is provided by the random operator r(p, t) : a(p, t) — > b(p, t) = a(p, t) + r(p, t). 
The operators (p) and (p) obey the following equations in 5?4 (see details in 
[3]): 

[w - K(p0 x \ P ) = F(p) + p(u> P ,e), (30) 
[w - X+(p)]6^ + (p) = F+(p) + p> P , e), (31) 

where p^ = (uj = p ,p)- Both equations (|30[) and (|3ip can be transformed into new 
equations for the frequency parts B^(x) and (x) 



idoB^ 



l(x)+ [ K(x-y)Bl l (y)dy = U(x), (32) 
id Bt + (x) + / - y) Bl + (y)dy = f+(x), (33) 



where 



W - / (27r) 3 2( ^+ w2) i/2 E4^)^W +P*(^P,e)]e^. (35) 



A=l 

3+ . 



The equations for field components B^{x) and B^ (x) (f3"2")) and (f3"3")) . respectively, are 
nonlocal within the presence of the formfactors K(x — y) and K + (x — y), respectively. 
In principle, these formfactors can admit the description of locality for nonlocal in- 
teractions. At this stage, it must be stressed that we have new generalized evolution 
equations (|32p and (|33|) . which retain the general features of the propagating and in- 
teracting of the quantum vector fields with mass m that are in the heat bath (thermal 
reservoir) and are chaotically distorted by other fields. For further analysis, let us 
rewrite the system of Eqs. ([3"2l and (|3"3"]) in the following form: 

id Bl(x) + K(x)*Bl(x) = f^x), (36) 
- l d B^{x)+K+{x)*Bf (x) =f+(x), (37) 
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where A(x) * B(x) is the convoluted function of the generalized functions A(x) and 
B(x). Applying the direct Fourier transformation to both sides of Eqs. (|36p and (|37|) 
with the following properties of the Fourier transformation 

F[K(x)*Bl(x)]=F[K(x)]F[Bj t (x)] (i = 1,2+), 

we get two equations 

\p° + k(p)]Bl(p)=F[U(x)], (38) 
{-p°-k + (p)]Bl + (p)=F[f+{x)}. (39) 
Finally, we have got the following equation for B^{p) field: 

[-p° + K+(p)][p a + K(p)}B^(p) = T», (40) 

where 

T» = {-P° + K+(p)]F[f M (x)] + [p° + K(p)]F[f+(x)]. 

We are now at the stage of the main strategy: one has to identify the field B^(x) 
and the random source operator Jr^ (x) , introduced in Eq. (|25[ with the Fourier 
transformed field B^ip) and T{p) in (jit})) , respectively. 

The next step is our requirement that Green's function G M „(p) in Eq. (|26|) and 
the function T^ijp), satisfying the equation 

{-p° + k + (p)}[p° + k(p)]f^(p) = (41) 

must be equal to each other, i.e. 

F[G ltu (p)-f ia ,(p)] = 0. 

The kernel operator K(p) is 

K{p) ~ 6^1 + ^, (42) 



where e = 2y kf + mf is the total energy of the final lepton-antilepton pair (with 

momentum ki and the mass mi for the lepton) produced within the decay of Z° 
boson being in the rest frame. To get k(p) in the form( |4"2"| we used the fact that 
the full Green's function G^ v {p) is given by the corresponding full Green's function 
of the scalar field [13] under the action by the differential operator {g^ u — m~ 2 p^p v ). 

5 Source size 

It has been emphasized [4] that there are two different scale parameters in the model 
considered here. One of them is the so-called "correlation radius" R introduced in ([5]) 
and ^ with (fT2"|) . In fact, this i?-parameter gives the pure size of the particle emission 
source without the external distortion and interaction coming from other fields. The 
other (scale) parameter is the stochastic scale L st which carries the dependence of the 
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particle mass, the a-coherence degree and what is very important — the temperature 
T-dependence: 



where 



1 



a(N)\p° - K(p)\ n(m,/3) 



4 



a{N) 4k? |1 - 5 k \ 2 n(m,l3) 



(43) 



and the lepton mass to; is neglected. 

It turns out that the scale L st defines the range of stochastic forces. This effect is 
given by a(iV)-coherence degree which can be estimated from the experiment within 
the two-particle BEC function C 2 (<3) when Q close to zero, C 2 (0), at fixed value of 
mean multiplicity (N) : 



a(N) 



C 2 (0) + v/2-C 2 (0) 



C 2 (0) = C 2 (0)/£(JV). 



C a (0) - 1 

In formula (|43[) . n(m,/3) is the thermal relativistic particle number density 



(44) 



(45) 



where iC 2 (...) is the modified Bessel function. 

The coherence function a is another very important one that summarizes our 
knowledge of other than space-time characteristics of the particle emission source, 
and the prediction of a from an experiment is very instructive aim itself. For a = 0, 
one actually finds 



KC 2 (Q)<£(A0(l + 7i 



-era- 



which is nothing other but the Goldhaber parameterization [11] with < 7 < 1 being 
a free parameter adjusting the observed value of C 2 (Q = 0). 

Within our aim to explore the correlation between Z°Z° the scale L st has the 
form 



12 a{N) kf (to 2 + M 2 ) (2IF 



2 (T_ \ 3 / 2 
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U-4 



(46) 



where the condition I (3 y/ to 2 + pt 2 > 1 for any integer I in (|45[l was taken into account. 
The only lower temperatures will drive L s t within formula (|46|) even if [i = and I = 1 
with the condition T < m. 

Note that the condition /1 < to is a general restriction in the relativistic " Bose-like 
gas" , and /i = to corresponds to the Bosc-Einstcin condensation. 
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For high enough Tno/i- dependence is found for L st : 



12C(3)a(7V)fcf T 3|l-4| 2 



(47) 




fj( m ,/3)^n(/3)~— C(3) 

7T Z 



. is taken into account. The origin of 



(48) 



where neither a Z° boson mass nor the fi - dependence occurred; £(3) = ' 3 = 

1.202 is the zeta- function with the argument 3. 

To be close to the experiment there is necessary to include transverse momenta, 
where the Z° boson mass m, in Eqs. P5|) . (U5J), (14"?| is replaced by the transverse 



Actually, the increasing of T leads to squeezing of the domain of stochastic force 
influence, and L st (T = To) = R at some effective temperature To. The higher tem- 
peratures, T > To, satisfy to more squeezing effect and at the critical temperature T c 
the scale L st (T = T c ) takes its minimal value. Obviously T c ~ 0(200 GeV) defines 
the phase transition where the chiral symmetry restoration will occur. Since in this 
phase all the masses tend to zero and a — ► at T > T c one should expect the sharp 
expansion of the region with L st (T > T c ) — ► oo. 

The qualitative relation between R and L st above mentioned is the only one we 
can emphasize in order to explain the mass dependence of the source size. 

6 Conclusions 

To summarize: the theoretical proposal for two-particle Bose-Einstein correlation 
function in case of Z° Z° pairs in pp collisions is carried out for the first time. 

The correlations of two bosons in 4-momentum space presented in this paper of- 
fer useful and instructive complimentary viewpoints to theoretical and experimental 
works in multiparticle femtoscopy and interferometry measurements at hadron collid- 
ers. 

We find the time dependence of correlation function calculated in time-dependent 
external field provided by the operator r(p,t) and the chaotic coherence function 
a(m,(3). The result can be compared with the static correlation functions (see, e.g., 
[14] and the references therein mainly devoted to heavy-ion collisions) and also can 
be used for experimental data fitting. 

The stochastic scale L st decreases with increasing temperatures slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is ap- 
proached. 

Our results first predicted for correlation radius R are bothZ boson mass and 
lepton energy dependent 



mass rriT — y^m 2 + pip. 



e m/5T 



(49) 



R 



- ,2/5 




m: 



3/10 J^ 3 / 10 
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for low values of To < m, while for higher temperatures, Tq > \J m 2 + /j, 2 , one has 



R ^ -• (50) 

a^\h\ Vb T^ 

The theoretical correlation radius R at temperature To decreases as Z°-boson mo- 
mentum increases. Both estimations (j49|) and (|50|) serve as the first approximation 
to explain the experimental data at different y/s and hence at T. We claim that the 
experimental measuring of R (in fm) can provide the precise estimation of the effec- 
tive temperature Tq which is the main thermal character in the Z°Z° pair emitter 
source (given by the effective dimension R) in the proper leptonic decaying channel 



Z°Z° — > llll with the final lepton energy y kf + mf at given a fixed by C^Q = 0) and 
(N) . Actually, Tq is the true temperature in the region of multiparticle production 
with dimension R = L st , because at this temperature it is exactly the creation of two 
particles (Z°Z°) occurred, and these particles obey the criterion of BEC. 
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